Neutrino Oscillation Effects in Indirect 
Detection of Dark Matter 



Nicolao Fornengo 

Dipartimento di Fisica Teorica, Universita di Torino 
and Istituto Nazionale di Fisica Nucleare 
via P. Giuria 1, F10125 Torino, Italy 



Abstract. If neutrino oscillation plays a role in explaining the atmospheric neu- 
trino deficit, then the same phenomenon would necessarily affect also the dark matter 
indirect-detection signal which consists in a muon-neutrino flux produced by neutralino 
annihilation in the Earth core. In this paper we investigate to which extent the upgoing- 
muon signal originated by neutralinos captured inside the Earth would be affected by 
the presence of i/^ — + Vr oscillation. 



1 Introduction 



Among the different techniques which have been proposed to search for dark 
matter particles detection of a neutrino flux by means of neutrino tele- 
scopes represents certainly an interesting tool. Although direct detection ||,^,^ 
at present appears to be somewhat more sensitive to neutralino dark matter 
nevertheless all the different possibilities are worth being explored. In 
this paper we discuss the flux of upgoing muons which are a consequence of i/^'s 
produced inside the Earth, with special emphasis on the role which is played in 
this kind of searches by the possible presence of neutrino oscillation. 



2 Upgoing /x's from neutralino annihilation in the Earth 

Neutralinos can be gravitationally captured inside astrophysical bodies Q , like 
the Earth and the Sun. Their subsequent annihilation can produce a flux of 
neutrinos, which then travel toward a detector located underground below the 
Earth surface. The differential flux, for each neutrino flavour i, is defined as 

where Fa denotes the neutralino annihilation rate, d is the distance between the 
detector and the source (which can be the center of the Earth or the Sun) , F is an 
index which lists all the possibilc final states which can be produced by neutralino 
pair-annihilation, B-^j denotes, for each final state F ^ the branching ratios into 
heavy quarks, r leptons and gluons. The differential spectra of neutrinos and 
antineutrinos generated by the r and by hadronization of quarks and gluons 
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and the subsequent semileptonic decays of the produced hadrons are denoted 
by dN ^_ JdE^. For more details, see for instance Refs. [p|P,p|,p^ . Here we only 

recall that the annihilation rate depends, through its relation with the capture 
rate of neutralinos in the Earth, on some astrophysical parameters, the most 
relevant of which is the local density pi . 

The best way of identifying the presence of these fluxes relies on the possi- 
bility to detect upward going muons inside a neutrino telescope. These upgoing 
muons would be produced by the component of the neutrino fluxes of Eq.(l). 
The charged-current interaction of the v^^^s with the rock below and close to 
the detector would produce a flux of muons. A double-differential muon flux is 
defined as 

where j = v^,v^, Na is the Avogadro's number, g{E^, E'^] X) is the survival 
probability that a muon of initial energy E'^ will have a final energy E^ after 
propagating along a distance X inside the rock and daj{E^, E'^)/dE'^ is the 
charged-current cross section for the production of a muon of energy E'^ from a 
neutrino (antineutrino) of energy i?^. 

A useful quantity for our discussion is the muon response function 

'""^ - dE,. , (3) 



dEu J dEfj^dEu 

where E^^ is threshold energy for detection of up-going muons. For Super 
Kamiokande and MACRO, E*''^ ~ 1.5 GeV Q. The muon response function 
identifies the neutrino energy range that is mostly responsile for the up-going 
muon signal. Fig. 1 shows a few examples of it, plotted as functions of the 
variable x — E,j/m^, where denotes the neutralino mass. Fig. 1 shows an 
approximate scaling of dN^/dE^, with m^. The maximum of the muon response 
occurs for neutrino energies of about Ei, ~ (0.4 — 0.6) m^, with a half width 
which extends from E'^ ~ 0.1 my- to E^ ~ 0.8 m^. 

The quantity which is actually measured is the total flux of up-going muons, 
which is deflned as 

/•™x fj AT 

can be calculated once a specific supersymmetric model is adopted. In the case 
of a model where all the supersymmetric parameters are defined and set at the 
electroweak scale (which we call here MSSM), the result for is shown in Fig. 2. 
We have varied the MSSM parameters in the ranges: 20 GeV < M2 < 1000 GeV, 
20 GeV < 1^1 < 1000 GeV, 90 GeV < ruA < 1000 GeV, 100 GeV < mo < 
1000 GeV, —3 < A < -1-3, 1 < tan/3 < 50. Up-to-date bounds and limits coming 
from accelerators and from BR(& — > 57) have been imposed. For a definition of 
supersymmetric models and their parameters, as well as the implementation of 
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the experimental limits on susy searches, see Ref. p\. For calculations of in 
supergravity inspired (SUGRA) models, see Ref. [^|. 

Fig. 2 also shows the present most stringent upper limit obtained by the 
MACRO Collaboration [|l3j . Super Kamiokande recently also reported a similar 
upper bound Q. 



3 Neutrino oscillation effect on the up— going muon signal 

The atmospheric neutrino deficit strongly points toward the indication that the 
i^fi may oscillate. The oscillation channel which best describes the anomaly is 

— J- vacuum oscillation [ p^|Jl^ . If this is the case, also the produced 
by neutralino annihilations would undergo an oscillation process. The range of 
energies involved in both atmospheric and neutralino-produced neutrinos is ap- 
proximately the same, while the baseline of oscillation of the two neutrino com- 
ponents is different Atmospheric neutrinos which induce upgoing muons cover 
a range of pathlengths which ranges from twice the Earth's radius, for vertical 
muons, to much shorter distances in the case of horizonthal muons. On the con- 
trary, neutrinos produced by neutralino annihilation in the central part of the 
Earth travel a distance of the order of the Earth's radius to reach the detec- 
tor. On the basis of the features of the I'f^ oscillation which are required to fit 
the experimental data on atmospheric neutrinos ||l5|,0 , we expect that also the 
neutrino flux from dark matter annihilation would be affected ||l7|] . 

For h'^ —>■ Vr oscillation, the i/^ flux is reduced because of oscillation, but 
we have to take into account that neutralino annihilation can also produce v-^ 
which in turn can oscillate into z/^ and contribute to the up-going muon flux. 
The muon neutrino flux can therefore be expressed as (we are considering only 
two-flavour oscillation) 

= p™^(y ^ + [1 - p™^(y ^ y )] , (5) 



where the vacuum survival probability is |18 



where Arr? is the mass square difference of the two neutrino mass eigenstates, 
is the mixing angle in vacuum and K is the Earth's radius. Fig. 3 shows the 
survival probability for two different values of the neutrino oscillation parameters 
which are inside their allowed ranges fl^Jl^ • Smaller (larger) values of Aw? have 
the effect of shifting the curves to the left (right). Comparing Fig. 1 with Fig. 3, 
we notice that the reduction of the up-going muon flux is stronger when there 
is matching between the the energy E]^ ~ 5.2 • 10^Z\m^(eV^) of the first (from 
the right) minimum of the survival probability and the energy Ey ~ O.Sm^^ 
which is responsible for most of the muon response in the detector. This implies 
that a maximum reduction of the signal could occur for neutralino masses of 
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the order of m^(GeV) ~ 10'^Z\m^(cV^). The Vr Vf_i oscillation makes the 
reduction of the muon flux less severe, but it is not able to completely balance 
the reduction effect because the original Vj- flux at the source is smaller than the 
flux. Therefore, the overall effect of the neutrino oscillation is to reduce the 
up-going muon signal. The upgoing muon flux for a neutralino in the MSSM, 
when neutrino oscillation is included, is given in Fig. 4. This, when compared 
with Fig. 2, shows the effect induced by the presence of oscillation. The ratio of 
the up-going muon signals in the presence and in the absence of oscillation is 
plotted in fig. 5. We notice that the strongest effect occurs for light neutralinos, 
since in this case the muon flux is mostly produced from neutrinos whose energy 
is in the range of maximal suppression for the oscillation phenomenon. The 
effect is between 0.5 and 0.8 for < 100 GeV. On the contrary, the fluxes 
for larger masses are less affected, and the reduction is less than about 20% for 
> 200 GeV. Figs. 2, 4 and 5 update the corresponding figures of Ref.||l^ by 
the inclusion of the new limits from accelerator on Higgs and Supersymmetry 
searches. 

4 Conclusions 

In this paper we have discussed to which extent neutrino oscillation can affect 
the up-going muon signal from neutralino annihilation in the Earth. While the 
experimental upper limit is, at present, practically not affected by neutrino oscil- 
lation , the theoretical predictions are reduced in the presence of oscillation. 
By adopting the neutrino oscillation parameters deduced from the fits on the 
atmospheric neutrino data ||l^,|6) , the effect is always larger for lighter neutrali- 
nos. For — > Vr the reduction is between 0.5 and 0.8 for < 100 GeV and 
less than about 20% for > 200 GeV. 
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Fig. 1. Muon response function dN^/dlogx vs. the parent-neutrino fractional energy 
X = E^/m^ for neutralino annihilation in the Earth. Different curves refer to different 
neutralino masses : = 50 GeV (solid), = 80 GoV (dotted), = 120 GeV 
(shot-dashed), = 200 GeV (long-dash), = 500 GeV (dot-dashed). 
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Fig. 2. Flux of up-going muons ^j^'^i'*'! from neutralino annihilation in the Earth, plot- 
ted as a function of m^. The solid line denotes the present upper limit from MACRO 
p^ ]. Crosses denote supersymmetric configurations for which the neutralino relic abun- 
dance fi-^h^ is larger than, or equal to, the value 0.05 (but not in excess of its cosmo- 
logical upper bound of 0.7 |)). Dots stand for Q^h^ < 0.05. The fp's are assumed not 
to oscillate. 
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Fig. 3. i/fi survival probability in the case of — > i/t oscillation. The solid line refers 
to sin^(26') = 1, the dashed line is for sin^(26') = 0.8. In both cases, Z\m^ = 5 ■ 10"^ eV 
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Fig. 4. Flux of up-going muons ^^'^'^^ ^ from neutralino annihilation in the Earth, plot- 
ted as a function of m^. The solid line denotes the present upper limit from MACRO 
p^ ]. Crosses denote supersymmetric configurations for which the neutralino relic abun- 
dance f2-^h^ is larger than, or equal to, the value 0.05 (but not in excess of its cosmo- 
logical upper bound of 0.7 |^]). Dots stand for H^h^ < 0.05. The f^'s are assumed to 
oscillate into i^r's, with oscillation parameters fixed at the best-fit values of Ref. fl^ : 
sin2(26i) = 1 and Am'^ = 3 • 10"^ eV^ 
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Fig. 5. Scatter plot of the ratio (^p''''*'')osc:iii/^J?'"'"' vs. tlie neutralino mass m^. 
(<p|^''^*'')osciii is tiie up-going muon flux in tlie case of Vfj, — > oscillation (shown 
in Fig. 4) , while (p^^'^^^ is the corresponding flux in the case of no oscillation and (plot- 
ted in Fig. 2). For the oscillation case, the neutrino parameters have been set at the 
best-fit values of Ref. sm^{2e) = 1 and Am^ = 3 ■ 10"^ eV 



